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Abstract: In this paper, we show that spin waves, the elementary excitation of the 
Heisenberg magnetic system, obey a kind of intermediate statistics with a finite maxi- 
mum occupation number n. We construct an operator reahzation for the intermediate 
statistics obeyed by magnons, the quantized spin waves, and then construct a corre- 
sponding intermediate-statistics reahzation for the angular momentum algebra in terms 
of the creation and annihilation operators of the magnons. In other words, instead of 
the Holstein-Primakoff representation, a bosonic representation subject to a constraint on 
the occupation number, we present an intermediate-statistics representation with no con- 
straints. In this realization, the maximum occupation number is naturally embodied in the 
commutation relation of creation and annihilation operators, while the Holstein-Primakoff 
representation is a bosonic operator relation with an additional putting-in-by-hand re- 
striction on the occupation number. We deduce the intermediate-statistics distribution 
function for magnons from the intermediate-statistics commutation relation of the creation 
and annihilation operators directly, which is a modified Bose-Einstein distribution. Based 
on these results, we calculate the dispersion relations for ferromagnetic and antiferromag- 
netic spin waves. The relations between the intermediate statistics that magnons obey 
and the other two important kinds of intermediate statistics, Haldane-Wu statistics and 
the fractional statistics of anyons, are discussed. We also compare the spectrum of the 
intermediate-statistics spin wave with the exact solution of the one-dimensional s = 1/2 
Heisenberg model, which is obtained by the Bethe ansatz method. For ferromagnets, 
we take the contributions from the interaction between magnons (the quartic contribu- 
tion), the next-to- nearest-neighbor interaction, and the dipolar interaction into account 
for comparison with the experiment. 
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1 Introduction 

At low temperatures, the elementary excitations of a magnetic system, a periodic system 
of localized spins coupled by exchange interaction, are magnons, the quantized spin waves. 
In this paper, we show that magnons obey a kind of intermediate statistics in which the 
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maximum number of particles in any quantum state is neither 1 nor oo, but equals a finite 
number n. That is to say, magnons are intermediate-statistics type quasiparticles. 

Let us first recall the common treatment of a magnetic system. Take the ferromagnet 
as an example. The Hamiltonian of the Heisenberg model of a ferromagnetic system 
describing the exchange interaction between neighboring spins reads 



H — ^ ^ Jmn^m ' Sn, (1) 



mn 



where Jmn is the exchange coefficient and Sm and S„ are the spins at m-th and n-th sites. 
For such a spin system, introducing the spin deviation operator at site i, 

N, = S-SI, (2) 

where 5| is the z-component of the spin operator and S| = 5" (5 + 1), we have 



S+ \N,) = ^/2S-iNe-l)VNe\Ne - 1) , 

S!\Ne)=S-Ne, (3) 

where \Ni) is the eigenstate of the spin deviation operator, i.e., Ni\Nn) = Nn\Ni), 
[5^,5^] = 2Sg5f^i, and [51,51^] = ±S^Si£' [1]. A natural restriction is 

0<Ne< 2S, (4) 

since Sf must be less than the magnitude of S^. The relation ([3]) leads to an operator 
realization for the angular momentum: 



S+ = ^2S - Neaf'"^ 
S^ = af"'^^y^2S^i, 

Sf = S- Ni, (5) 
where af°'^'^, af°''^'^\ and N^°^^ satisfy the bosonic commutation relation: 

Base Bose\^ _ r 

g , a^, — o^£' , 

'Ne,afr'^'\ = af'^H,^ and [A^£,af°^^] = -af°^'=5„,. (6) 

Prom equation ([6]), we can see that a^°^^'^ creates and a^°^^ annihilates a localized spin 
deviation at a single site. This is the Holstein-Primakoff representation of angular mo- 
mentum algebra [2]. 

It is known that the Holstein-Primakoff representation is not a genuine bosonic re- 
alization of angular momentum algebra, since though the operators satisfy the bosonic 
commutation relation ([6l), there still exists an additional restriction on the value of the 
spin deviation equation When > 25, the representation is not faithful, while 
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in the Bose-Einstein case, Ni can take on any value. That is to say, the Holstein-Primakoff 
representation corresponds essentially to a kind of intermediate statistics with a maximum 
occupation number 25". In the Holstein-Primakoff representation, though there exists a 
maximum occupation number n = 25, the maximum occupation number is not embodied 
in the operator relation ([6]); in fact, the Holstein-Primakoff representation is a bosonic 
realization with a putting-in-by-hand restriction on the occupation number. As a result, 
however, when using the Holstein-Primakoff representation to solve the spectrum, only 
the operator relation is used, so the influence of the restriction on the occupation number 
is ignored. 

The above analysis shows that the spin deviation does not obey Bose-Einstein statistics, 
but obeys a kind of intermediate statistics with a finite maximum occupation number. In 
fact, as discussed in [3lll], if one wants to construct an operator realization for the angular 
momentum algebra by a single set of creation and annihilation operators, he needs a kind 
of intermediate statistics rather than the Bose-Einstein or Fermi-Dirac case. Nevertheless, 
in the common treatment of spin waves, the spin deviations are regarded as bosons: the 
commutation relation of creation and annihilation operators is taken as equation ([6]) and 
the statistical distribution is taken as the Bose-Einstein distribution. 

In this paper, we will construct an intermediate-statistics operator realization for the 
angular momentum algebra, in which the maximum occupation number is naturally em- 
bodied in the commutation relation of creation and annihilation operators. Therefore, all 
results based on this intermediate-statistics realization can naturally take the influence 
of the restriction of occupation number into account. As a comparison, the Holstein- 
Primakoff representation is a constrained bosonic representation, i.e., a bosonic realiza- 
tion subject to a constraint on the occupation number, so the influence of the maximum 
occupation number cannot be taken into account directly. 

From the intermediate-statistics realization for the angular momentum algebra, we 
calculate the corresponding intermediate-statistics distribution function. We show that 
the statistical distribution that magnons obey is a modified Bose-Einstein distribution. 

Based on the intermediate-statistics realization for the angular momentum algebra and 
the intermediate-statistics distribution function, we calculate the dispersion relations of 
spin waves for ferromagnet and antiferromagnet. 

The magnons, as shown in the present paper, obey a kind of intermediate statis- 
tics. As comparisons, we will discuss the relations between the intermediate statistics 
that magnons obey and the two important kinds of intermediate statistics: Haldane-Wu 
fractional statistics [5] and the fractional statistics of anyons [6j. Haldane-Wu fractional 
statistics is constructed based on the generalization of the Pauli exclusion principle. The 
concept of anyons is introduced by analyzing the symmetry properties of the wavefunc- 
tion of identical particles: the change of the phase factor of the wavefunction when two 
identical particles exchange, instead of or —1, is generalized to an arbitrary phase 
factor e*^. Each of these two kinds of intermediate statistics has an intermediate-statistics 



4 



parameter: g (Haldane-Wu) and a = 9/tt (anyon). The roles of g and a in Haldane-Wu 
statistics and in the fractional statistics of anyons are just as the role of the parameter 
n in the intermediate statistics obeyed by magnons. In this paper, we will discuss the 
relations of g and n and a and n, respectively. 

A special case of the Heisenberg model, the one-dimensional s = 1/2 Heisenberg spin 
chain, can be solved exactly by the Bethe ansatz [7l[8l[9]. By this exact solution, we can 
check the validity of our result directly. Concretely, we will compare the exact spectrum 
obtained by the Bethe ansatz method with the spectrum obtained by the intermediate- 
statistics method and the spectrum obtained by the Holstein-Primakoff method. 

Moreover, we also compare our result with the experimental data of EuO given in |10j . 
The result shows that at low temperatures and low frequencies, the intermediate-statistics 
spin wave model is more accurate than the bosonic spin wave model. 

The present paper discusses the intermediate statistics and the intermediate-statistics 
realization of angular momentum algebra and their applications to magnetic systems. On 
the one hand, the realization of angular momentum algebra has been of interest for a long 
time. Besides the well-known Schwinger representation [11] and the Holstein-Primakoff 
representation [2], there are many other schemes, including the realization of su(2) alge- 
bra O m [121 131 E]) and the realization of suq (2) and suq (n) algebra [151 US]- On the 
other hand, the spin wave plays an important role in magnetic problems \17\ I18j . The 
concept of spin wave has become a widely applied tool in the fields related to magnetism. 
It has been applied to study magnetic semiconductors [TO], quasiequilibrium spin systems 
[20j . ballistic thermal transport [21] and the thermodynamics (22j in the Heisenberg spin 
chain, one-dimensional ferromagnetic Bose gases [23], spin- wave excitations in cylindrical 
ferromagnetic nanotubes [Ml? etc. In experiment, the properties of spin wave have been 
directly measured [251 126j. The realizations of angular momentum algebra are successful 
in describing magnetism in various quantum systems [27J. In the application of the real- 
izations of angular momentum algebra to magnetism, the realization is either bosonic or 
fermionic. However, as shown above, magnons obey neither Bose-Einstein nor Fermi-Dirac 
statistics. Therefore, a kind of intermediate-statistics treatment is needed. As generaliza- 
tions of Bose-Einstein and Fermi-Dirac statistics, many schemes of intermediate statistics 
have been discussed [3 [6l [281 [291 [30] . 

In section [2] we construct an intermediate-statistics operator realization for the angu- 
lar momentum algebra. In section [S] we calculate the intermediate-statistics distribution 
function for magnons based on the commutation relation between creation and annihila- 
tion operators of magnons. In sections [3| and [5l we calculate the dispersion relations for 
ferromagnetic and antiferromagnetic spin waves. In section [Ul we discuss the relation be- 
tween the intermediate statistics that magnons obey and Haldane-Wu fractional statistics 
and the relation between intermediate statistics and the fractional statistics of anyons. In 
section [71 we compare our result with the exact solution obtained by the Bethe ansatz 
method. In section [51 we compare the dispersion relation of a ferromagnetic system, which 
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is calculated based on intermediate statistics, with the experimental data of EuO. The 
conclusions and discussions are given in section [H 



2 Intermediate-statistics operator realization for angular mo- 
mentum algebra 



For solving the spectrum of a magnetic system with Hamiltonian ([T]), we need a rep- 
resentation of angular momentum algebra. As discussed above, the Holstein-Primakoff 
representation is a bosonic realization with a putting-in-by-hand restriction on the occu- 
pation number, and when using it to solve the spectrum, since only the operator relation 
is taken into account, the information of the restriction on the occupation number is ig- 
nored. For taking the influence of the restriction on the occupation number into account, 
we need a representation in which the information of the maximum occupation number 
is embodied in the operator relation rather than put in an additional restriction by hand. 
In this section, we construct an intermediate-statistics operator realization for the angular 
momentum algebra, in which the maximum occupation number is naturally embodied in 
the commutation relation of creation and annihilation operators. 

In the case of the spin wave, we only focus on the low-lying excitation. In other words, 
in our case the spin deviation is always very small, and then equation ^ can be 
expanded around = 0; only taking the next-to-leading-order contribution into account, 
we have 



St \Np 



S7 \Nf, 



S'i \No 



2S 1 



1 



AS 



^,\N,-l) 



V2SVNe + 1 (l - ^) |A^^ + 1) 



S-Nf. 



(7) 



For constructing a realization of angular momentum algebra, we first introduce an 
intermediate-statistics operator realization: 



ae,al, 



1 



Ne 
n 



1 _ El 

^ 2n 



1 



Nt 
n 



1 



n 



2n J 



Ne-1 
2n 



N£,a\i = a^dfji and [N^^a^i] = —afd^^', 



(8) 



where n = 25. It can be directly checked that such an operator realization corresponds 
to a kind of intermediate statistics with a maximum occupation number n. 

By creation, annihilation, and number operators, a£, o!, and A'^, we can construct an 
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intermediate-statistics realization of angular momentum algebra: 

<?+ 



Si 



S -Ne. 



It can be directly checked that, with the commutation relation 
the operator relation of angular momentum: 



(9) 

) , , , and 5| satisfy 



(10) 



[sf,sf] 

[5/, S"^,] 

In this scheme, a\ and ai are the creation and annihilation operators of a localized spin 
deviation at a single site, and Ni is the spin deviation operator which can be expressed as 



a\ai + {n + l) - J {n + lf' - (2n + 1) a\ai 



(11) 



a^a£ + 4n 

It should be emphasized that in such intermediate statistics, ^ a\ai. Concretely, for 
a given set of creation and annihilation operators, a\ and a^, the corresponding number 
operator Ng^ can be constructed from the operator relations of a\, a\, and A'^. As shown in 
[3] and , the number operator of intermediate statistics in general cannot take the form 
of A'^ = a\ai. The realization of the number operator (llip shows that the spin deviation 
corresponds to a kind of intermediate statistics. 

For describing the nonlocalized excitations of such a magnetic system, taking transla- 
tional symmetry into account, we replace the creation (annihilation) operator a\ {an) which 
creates (annihilates) localized spin deviations with the creation (annihilation) operator 
(5k) which creates (annihilates) nonlocalized excitations by the transformation 

1 



E 

k 

E 



where W is the number of lattice sites. Then 



1 



n 



1 _ ilk 

2n 



+ iVk 



1 



e^^-^6k, 



-ik.£,t 
e "k' 



n 



(12) 



1 _ 

^ 2n 



2n 



3kk'> 



^k'^kk' and [iVk,^k'] = -f'k4k' 



(13) 



-^k) and 6k are the number, creation, and annihilation operators of magnons, respec- 
tively, which describe the nonlocalized elementary excitation of the system. 
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The above result shows that magnons, the elementary excitation of a magnetic system, 
obey the intermediate statistics defined by (|13p . which has a maximum occupation number 
n. Only when n oo, such intermediate statistics returns to Bose-Einstein statistics. 
In other words, the spin wave is essentially an intermediate-statistics type elementary 
excitation. 



3 The intermediate- statistics distribution function 



When calculating the spectrum of a magnetic system, one needs to use the statistical 
distribution function of the magnon. In the common treatment, the statistical distribution 
is approximately taken as the Bose-Einstein distribution. However, as discussed above, 
the magnon in fact obeys intermediate statistics. In this section, we seek for the statistical 
distribution for an ideal magnon gas based on the commutation relation (I13p . 
The average particle number of k state can be obtained by 



1 

z 



-.Tr 



(14) 



where Z = Tre~^^^~'^^\ N is the total particle number operator, and is the chemical 
potential. For the low-lying excitation, we can expand the expression of the particle 
number operator, similar to equation (jlip . as 



(2n + l) 



{2n+lf 



An {n + 1) "^""^ 16n2 (n + 1)^ V^l^^^ 



By the relation 



from (I14D. we achieve 



(15) 
(16) 



Z An{n+ 1) 



X <^ Tr 



2n + l 

4n {n + ly 



tTt 



(17) 



where ek is the energy of k state. Based on the operator relation obtained in the above 
section, we construct 



1 



n 



2n 



2 and bvbl = (iVk + 1) 



1 



n 



2n ) 



(18) 



Ignoring {N^), we can solve (A'^k): 

im = - 



1 _ 4n^+4n^+2n+l 
16n2(n+l)^ 



/3(ek-M) 



1 



6n3+8ra^+4n+l 
8?i2(n+l)^ 



(19) 
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It can be directly seen from this equation that the statistical distribution defined by the 
commutation relation (113p is a modified Bose-Einstein distribution, and when n — > oo, 
equation (|19p returns to the Bose-Einstein distribution. In other words, the spin wave 
obeys a modified Bose-Einstein distribution. 



4 Intermediate-statistics ferromagnetic spin waves 

In this section, we calculate the dispersion relation of a ferromagnetic spin wave. The 
Hamiltonian of the Heisenberg model reads 



H = — Ji ^ • — J2 ^ • H 



(20) 



where 5i and 82 connect spin £ with its nearest and next-to-nearest neighbors and Ji and 
J2 denote the exchange parameters corresponding to the nearest-neighbor and next-to- 
nearest-neighbor couplings. 

We first consider the nearest-neighbor contribution. The Hamiltonian reads 

1 



(21) 



Substituting the intermediate-statistics representation of angular momentum algebra ([9]) 
into (j2ip . we have 

H = Ho + H2 + H4 (22) 



with 



1,6^ 

n 



£,Si 



(23) 
(24) 



and 



(25) 

where Zi is the number of nearest neighbors. 

In the following, by the above operator relations and the statistical distribution func- 
tion given in sections [2] and [3l we can calculate the spectrum directly. 



4.1 The quadratic contribution 

First, we calculate the contribution from the terms quadratic in the creation and annihi- 
lation operators. Substituting (fT2|) into (p^ gives 



H2 = JiZm ^ iVk - JiZi^ Yl Ubu + 



(26) 
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Here 7^ is defined to be 7^ = Ylsi e^^'^^ with 7^ = ^^^^ due to the symmetry. Then 
substituting (fTSj) into (f26]) and preserving only the first-order contribution of A^k gives 

1 



H2 = JjZiny^ 

k 



2 (2n + 1 



^7k^ 



iVk. 



(27) 



We achieve the dispersion relation of magnons to the second order: 



(2) 



JiZin 



1 



..-7k^ 



(28) 



2(2n + 1) 

The second term in ()28p is the modification coming from the influence of the restriction 
on the occupation number. 

4.2 The quartic contribution: the interaction between magnons 

The contribution from the terms quartic in the creation and annihilation operators de- 
scribes the interaction between magnons. A similar treatment can also be used to deal 
with the quartic terms. 

From equation (j25p . by the operator relations given above, up to quartic terms, we 
obtain 

2 



Ha 



-'4 E 



w 
1 + 



kik2k3k4 
1 



1 + 



1 



(^ki -k2+k3 ,k4 7k^3-k4 ^ki ^k2 ^k* 



An {n + 1) 

'5ki-k2-k3,-k4 7^1^ki^'k2&k3?'k4 + 5ki-k2-k3,-k4 7kl^ki&k2&k3^k4 

(29) 



4n (n + 1) 

+<5ki+k2-k3,k4 7ki^ki^k2^k3&k4 + 4i+k2-k3,k4 7kl^ki^k2^k3^k4) } 

For long- wavelength spin waves, the main contribution comes from the interactions that 
do not change the state of the spin wave. For example, for the term in proportion to 
^ki^k2^k3&k4) only the contributions corresponding to ki = ka and k2 = k4, or ki = k4 
and k2 = ka remain [9j. Only taking these contributions into account, we approximately 
achieve 

w 



kik2 



loblM^lb^^ + 7k2-k,&Lfek4 6k2bL 



1 + 



1 



4n (n + 1) 



An (n + 1) 

ll^AAAb^^ + 37^2&lc,ftki6k2&L + (37^1 + 7k2) blb^A.bi., 



By equation (fTSj) . up to N^, we obtain 
Z 



^ = -^4^[(^-<Hl-7k"2) 

kik2 

7k^,+7k^,+27k^,7k^ 



+ 



2(2n + 1) 



+ 



7k^.7k^, 



An {n + 1) (2n + 1) 



iVkiiVk2- 



(30) 

(31) 
(32) 
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In the calculation, the relations 7^ = 1, = TkaTki [3ll> Ek 7k = 0' Ekika ^ki^^i 

Ekika lk2^^i^i^2 have been used. 

When the number of excited magnons fluctuates little, we can take the approximation 

m 



(33) 



Then 



-2^4 E 



'1 -.A-) + 2i±ijL^« 



+ 



qk 

7q^7k^ 



+ 



4n (n + 1) (2n+ l)' 
7^,+7^, + 27K 



kik2 



7k".7k", 



2(2n + 1) 



4n (n + 1) (2n + 1) 



(34) 



Consequently, the contribution from the interaction between magnons to the dispersion 
relation reads 



(4) 




7^ 



' 2(2n + l)2 



7^ 



2(2n + l)^ 4n(n + l) 



(35) 



It should be emphasized that the statistical distribution function (A'^q) is the intermediate- 
statistical distribution given by , rather than the Bose- Einstein distribution as that in 
the common treatment. 



4.3 Comparison with the result of the Holstein-PrimakofF representation 

From the above discussion on the commutation relation of creation and annihilation op- 
erators, we have already known that magnons obey a kind of intermediate statistics with 
a maximum occupation number n, and when n — > 00 this intermediate statistics returns 
to Bose-Einstein statistics. In the common treatment, one approximately assumes that 
magnons obey Bose-Einstein statistics. Under this approximation, the dispersion relation 
reads [9] 



hio^ = J,Z,2S (1 - 7^) - 2Ji|i (1 - 7^) 



q q 



(36) 



where (^N^"^'^'^ denotes the Bose-Einstein distribution. 

Comparing the dispersion relation (|36p with the dispersion relation ()28p and (|35p . 
we can see that when replacing Bose-Einstein statistics by intermediate statistics, some 
additional modifications relying on the maximum occupation number n appear. When 
the maximum occupation number n — > cxd, such modifications vanish. 
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Especially, the statistical distribution function (^N^"''^'^'j in the Holstein-Primakoff re- 
sult ()36p is the Bose-Einstein distribution, while in the present result (|35p . the statistical 
distribution function (A'^q) is the intermediate-statistics distribution (|19p which is a mod- 
ified Bose-Einstein distribution. 



4.4 Other contributions 

For comparison with the experiment, we need to consider all effects as possible. In this 
subsection, we discuss the contribution from next-to-nearest neighbors and dipolar inter- 
actions. 

4.4.1 The next-to- nearest-neighbor contribution 

The contribution from the next-to-nearest-neighbor coupling can be calculated directly by 
the same procedure: 




2{2n + iy 4n(n + l) 



where Z2 is the number of next-to-nearest neighbors and = ^^^^^ ■ 



(37) 



4.4.2 The dipolar interaction 

In the above, we only take the contribution from the exchange coupling into account. 
Besides the exchange coupling, there still exists a classical dipolar interaction which is 
caused by the interaction between the magnetic moments. The Hamiltonian of the dipolar 
interaction reads 



1 



3 {Si 



(S, 



(38) 



where is the Bohr magnon and g = 2 is the Lande factor. In principle, we need to 
substitute the intermediate-statistics operator realization of angular momentum algebra 
([9]) and ([8]) into this Hamiltonian, and, then, calculate the influence of the intermediate- 
statistics dipolar interaction to the dispersion relation. However, since the contribution 
from the dipolar interaction, in comparison with the contribution from the exchange cou- 
pling, is small, we ignore the intermediate-statistics modification to the dipolar contri- 
bution, i.e., for the dipolar interaction, we approximately use the result obtained by the 
Holstein-Primakoff representation. 
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Under the assumption of isotropy, the contribution from the dipolar interaction, based 
on the result given in [32], can be approximately expressed as 

^k^ = f 5/"bM(T), (39) 
where M (T) is the magnetization. Then the dispersion relation reads 

^0;^°*'^' =nuj^ + hw^'P, (40) 
where huj]^ comes from the exchange coupling and hw'^^ comes from the dipolar interaction. 



5 Intermediate-statistics antiferromagnetic spin waves 



It has been shown that there do exist quantized spin waves in antiferromagnets [18]. In 
this section, we calculate the dispersion relation for antiferromagnetic spin waves based 
on the intermediate-statistics scheme. 

For antiferromagnets, the spin structure of the crystal is considered as two interpen- 
etrating sublattices A and B with the property that all nearest neighbors of a spin on A 
lie on B, and vice versa. The Hamiltonian reads 



H = 2JY^ 



es 



(41) 



where £ runs over all sites of sublattice A. In this paper, we only consider the contribution 
from the nearest neighbors. 

The antiferromagnetic ground state is approximately taken as the Neel state, in which 
the z component of each spin is S in sublattice j4, and —S in sublattice B. 

The excited state of antiferromagnets can be treated by the similar treatment of ferro- 
magnets. Similar to equation ([9|), introduce two sets of operator realizations for sublattice 
A and B, respectively: 

S^^ = V2Sa\ (l - ^] 



AS 



S%, = S-Ne, (42) 

and 

5+^ = V2Scl f 1 - ^ 



4S 



S%, = -{S-N,), (43) 
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where aj, and c^, q are the creation and annihilation operators of the spin deviations 
on sublattice A and satisfying the operator relations ([8]) and ([TT]) . respectively. 
Substituting (^2]) and (03]) into Hamiltonian ([IT]) gives 



= 2 J ^ {(5 - A^£) (-5 + N,+5 

Ne+s 



+S 



45 



45 



45 y 



(44) 



For low-lying excitations, we only take the contribution from the terms quadratic in the 
creation and annihilation operators into account: 



Here W denotes the number of the sites of sublattice A. 
Introduce the transformations 



es 



H ~ -2JWZS'^ + 2 J5 {Ne + Ne+s) + 2 J5 ^ (aece+s + alc^ 



at 



-y 



e *^'^6k and a! 



-y 



(45) 



(46) 



and 



dk and c 



(47) 



Using the relation between Ni and a\, ai, cj, and q and substituting the above transfor- 
mations into (|45p. up to the quadratic contribution, gives 

H -2JWZS' + ?^^X) JSZ (blbu + + 2JSZ ^ 7k (&kd-k + bldl^) . 

^ k k 

(48) 

For diagonalizing the Hamiltonian (|48|) . we introduce the Bogoliubov transformation 
which mixes the operators of the two sublattices: 



&k = -"kOk + Wk/3k, c?-k = iik/3k + ^^kal, 
ftt = Ukaj, + Uk/?k, c^lk = ""k/^k + ^^kOk- 



(49) 



It can be checked directly that the Hamiltonian (|48p can be diagonalized when Uk and Vk 
are taken as 



2 1 
k 2 » 



- < 



16n2 (n + 1)2__ 2 
(2n + 1) 



7k 



-1/2 



+ 1 



16n2(n+l)^ 2 
(2n + l)^ 



-1/2 



(50) 
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Then, from equation (I48p . ignoring the high-order contribution, we achieve 

(2n + l)2 



H = -2JWZS^ + 2JSZ 



+ 



4n (n + 1) 
By the operator relation 



2n + lY 2 



in (n + 1) 

ttkOk + /5k/?k 



■"k + Tk^k^k 



"kOk + /3k/?k 



(51) 



1 



2n 



and ak^k ~ i^ak + 1) 



?1 



1 



Nek 
2n 



(52) 



(the operator relation for (3k is the same as that of Ok), ignoring the high-order contribu- 
tion, from (1511). we achieve 



H = -2JWZS'^ + iJSZ 



{2n + lf 2 



4n (n + 1) 



% + Tk^^k^^k 



+ ^/ia;k(iVak + iV/3k), (53) 



where the dispersion relation for antiferromagnetic magnons is 

1 . 1 



hu!k = 2JSZ 



Uk + vi + 27kUkVk + 



-2'l]iUkVk 



(54) 



An{n + l) (2n + l) 

Similar to the ferromagnetic case, the antiferromagnet magnon obeys intermediate 
statistics rather than Bose-Einstein statistics. When approximately regarding the magnon 
as bosons, based on the Holstein-Primakoff representation, the dispersion relation of the 
antiferromagnet magnon reads 



huj 



HP 



(55) 



where the superscript "HP" denotes that the corresponding result comes from the method 
of the Holstein-Primakoff representation. Comparing the dispersion relation (|54p with 
(|55p . we can see that when the maximum occupation number n — > oo, our result returns 
to the result of the Holstein-Primakoff representation which regards magnons as bosons. 
The magnitude of the modification relies on the value of n. 



6 Comparison with other schemes of intermediate statistics 

6.1 Comparison with Haldane-Wu fractional statistics 

The above result shows that the magnons obey a kind of intermediate statistics with 
the statistical distribution (|19|) . Let us compare this statistical distribution with another 

kind of intermediate statistics Haldane-Wu fractional statistics [5j. The Haldane-Wu 

distribution function reads 
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where uj is determined by g In (1 — w) — Ino; = /3 {e^ — fi). 

In the Haldane-Wu distribution, there is an intermediate-statistics parameter g, and 
in the intermediate statistics obeyed by magnons, the intermediate-statistics parameter 
is n. A relation between these two intermediate-statistics parameters can be obtained by 
comparing the second virial coefficients. The second virial coefficient of a z^-dimensional 
ideal magnon gas with the dispersion relation e ck can be obtained directly: 

1 i^Tj^) 4n3 + 8n^ - 2 
02 2T (f + 1) 4n3 + 8n2 + 2n - 1 ■ ^ ^ 

The second virial coefficient of an ideal gas obeying Haldane-Wu fractional statistics reads 
(^2^ = — (1 — 2g) /2'^/*+^ [33j. Comparing these two second virial coefficients gives 



i^r(f ) in-^ + 8n^ - 2 
2T + l) 4n3 + 8n2 + 2n - 1 



(58) 



6.2 Comparison with the fractional statistics of anyons 



It is also interesting to compare this intermediate statistics with the fractional statistics of 
anyons, another scheme of intermediate statistics [6]. For the case of anyon, we of course 
only focus on two dimensions. 

The concept of anyons is introduced by generalizing the change of the phase factor 
of a wavefunction when two identical particles exchange to an arbitrary phase factor e*^. 
6 = and 9 = ir correspond to Bose-Einstein and Fermi-Dirac cases, respectively. 

The second virial coefficient of an anyon gas reads 



02 = (1 -4a + 2a2) , (59) 

where a = O/ir. Comparing this result with the second virial coefficient ()57p with v = 2 
and s = 2 gives 

a = 1 - Jl- , '^"'^-^ (60) 
V 8n3 + 16n2 + 4n - 2 ^ ^ 

7 Comparison with the exact result of the Bethe ansatz 
method: the spectrum 

In this section, we compare our result with the exact solution of the one-dimensional spin 
1/2 Heisenberg model. 

By the Bethe ansatz, one can find the exact solutions of certain one-dimensional quan- 
tum many-body models. Taking ferromagnets as an example, we compare our result given 
in section m with the exactly solved one-dimensional s = 1/2 Heisenberg model with two 
down spins. 
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I ■ 





I s 



— I 1"^ 

□ Holstein-Primakoff 

* intermediate statistics 

• Betlie ansatz 




ka 



Figure 1: The spectra of the one-dimensinal s = 1/2 Heisenberg chain given by the 
Bethe ansatz method (the exact solution), by the Holstein-PrimakofF method, and by the 
intermediate-statistics method. In the unshadowed area, in comparison with the exact 
resuh, the intermediate-statistics result is better than the result given by the Holstein- 
PrimakofF method. 
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The exact spectrum of the one-dimensional s = 1/2 Heisenberg model with two down 
spins is given by [71 [U [9] 



E — Eq = J {2 — cos kia — cos k2a) 



(61) 



where a is the lattice constant and ki, k2 are determined by 



Nkia 



27rAi + 9 




(62) 



with Ai, A2 = 0,1,2, - • • , A^ - 1 and A2 > Ai. 

Moreover, our result of the spectrum for the corresponding case can be directly ob- 
tained by equation (p8]) with n = 1. 

The spectra obtained by the Bethe ansatz (the exact one), by the Holstein-Primakoff 
method, and by the intermediate-statistics method are sketched in figure 1. 

Comparing with the exact result obtained by the Bethe ansatz, we can see that in 
some cases (the unshadowed area) our result (the intermediate-statistics magnons) is more 
accurate than the standard Holstein-Primakoff result (the bosonic magnons). 

8 Comparing with the experiment 

From equation (|4U|) . we can obtain the relation between the spin- wave energies and the 
temperature by the self-consistent calculation. We will consider the spin-wave dispersion 
relation of EuO since the spin-wave dispersion for EuO is isotropic. The Ev?^ ions in 
EuO form simple fee lattices, so the number of the nearest neighbors and next-to-nearest 
neighbors are Zi = 12 and Z2 = 6, the exchange parameters Ji to nearest neighbors and 
J2 to next-to-nearest neighbors are Ji = 0.606A;b and J2 = O.llQfes [35], where ks is the 
Boltzmann constant, and S = 7/2. 

The calculation results are plotted in figure 2. The experimental data are taken from 



In comparison with the experimental data, we can see that at low temperatures and 
low frequencies, the result of the intermediate-statistics spin waves is more accurate than 
the result of the bosonic spin waves, and at high temperatures and high frequencies, the 
result of bosonic spin waves is better. 

9 Discussion and Conclusions 

It is shown that magnons, the elementary excitation of a Heisenberg magnetic system, 
obey a kind of intermediate statistics with a maximum occupation number n = 2S. In 
the common treatment, the solution of the spectrum of a magnetic system is based on the 



m- 
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10 20 30 40 50 60 70 

T(K) 



Figure 2: The spin- wave energies in EuO. The dotted hnes represent the intermediate- 
statistics result and the sohd hnes represent the result of the Holstein-Primakoff represen- 
tation. The experimental data are taken from [lOj . 
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Holstein-Primakoff representation which is a bosonic operator relation with an additional 
restriction on the occupation number. Since the information of the maximum occupation 
number is not embodied in the operator relation, the influence of the restriction on the oc- 
cupation number is not reflected in the result of the spectrum. Consequently, the magnons 
are approximately treated as bosons in the Holstein-Primakoff treatment: the commuta- 
tion relation of creation and annihilation operators is taken as the bosonic commutation 
relation and the statistical distribution is taken as the Bose-Einstein distribution. 

In this paper, we construct an intermediate-statistics operator realization in which 
the information of the maximum occupation number which is equal to an integer n is 
embodied in the commutation relation of creation and annihilation operators rather than 
putting in a restriction on the occupation number by hand. Then, from the operator 
relations, we directly deduce the corresponding statistical distribution function, which 
is a modified Bose-Einstein statistical distribution and will return to the Bose-Einstein 
distribution when taking the maximum occupation number n to be oo. 

It is the starting point that there is a natural relation between the angular momen- 
tum and the intermediate statistics with a given maximum occupation number. For the 
intermediate statistics with a maximum occupation number n, there are n + 1 states, |0), 

|2),- • • , \n). For the angular momentum S, there are 25 + 1 states |— 5), \—S + 1),- • • , 
\S — 1), \S). This naturally leads us to relate the n + 1 states, |0), • • , |n), to the 
25 + 1 angular momentum states, |— 5), |5). Consequently, we have the relation 
n + 1 = 25 + 1, and then n = 25. From this, we can construct an intermediate-statistics 
realization and reveal that the statistics of magnons is intermediate statistics. 

Based on the results of the intermediate statistics provided in sections [2] and [3l we cal- 
culate the dispersion relation of the ferromagnetic spin wave up to the quartic contribution, 
in which the influence of the interaction between magnons is taken into account, and the 
dispersion relation of the antiferromagnetic spin wave up to the quadratic contribution. 
Compared to the result of the Holstein-Primakoff representation, the bosonic operator 
relation is replaced by the intermediate-statistics operator relation, and the Bose-Einstein 
distribution is replaced by the intermediate-statistics distribution, so the influence of the 
restriction on the occupation number is naturally taken into account. Moreover, we also 
take into account the next-to-nearest-neighbor contribution and the influence of the clas- 
sical dipolar interaction which is caused by the interaction between the magnetic moments 
in the ferromagnetic case. 

Magnons obey a kind of intermediate statistics. As comparisons, we discuss the rela- 
tions among the intermediate statistics obeyed by magnons, Haldane-Wu fractional statis- 
tics, and the fractional statistics of anyons. The relations among the three intermediate- 
statistics parameters are given. 

For discussing the validity of our result, we compare our result with the exact solution 
of the one-dimensional spin 1/2 Heisenberg model obtained by the Bethe-ansatz method. 

Our results of the dispersion relation of the magnetic systems are based on intermediate 



20 



statistics, in which the maximum occupation number is an integer n equahng 25. We 
compare our result with the result by the Holstein-Primakoff representation in which 
magnons are assumed to obey Bose-Einstein statistics and with the experimental data of 
EuO. The result compares well with the experiment. 

In a word, the elementary excitation of the Heisenberg magnetic system obeys a kind 
of intermediate statistics with an finite maximum occupation number n = 2S rather than 
Bose-Einstein statistics. 
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